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I.  INTRODUCTION  AND  SUMMARY 


In  ballistics  work,  and  especially  in  small  arms  firings 
at  targets  to  determine  accuracy,  it  is  quite  natural  to 
measure  the  closeness  of  a group  of  shots  by  the  "extreme  spread"; 
or  the  greatest  distance  between  any  two  shots  of  a group. 
Moreover,  the  extreme  spread  may  be  determined  quickly  with 
a ruler  and  does  not  require  any  detailed  or  involved 
computation,  as  does  the  round- to -round  standard  deviation 
in  each  direction,  or  the  mean  radius  of  the  shots  of  a 
group,  for  example.  It  is  for  these  reasons  that  ballisticians, 
riflemen  and  others  have  long  had  a great  interest  in  the 
extreme  spread,  for  it  is  truly  the  most  rapid  measure  of 
dispersion  of  shots  on  a target.  Sometimes  the  extreme 
spread,  or  the  maximum  distance  between  pairs  of  points  on 
the  target,  is  called  the  "group  diameter",  but  there  is  a 
very  subtle  difference  between  the  two  when  one  delves  into 
the  general  problem  on  a statistical  basis.  We  do  not 
intend  to  cover  all  the  pertinent  details  relating  to 
the  statistical  analysis  of  patterns  of  shots  on  a target  here, 
but  interested  readers  might  well  study  the  booklet  of 
Grubbs  [1964].  Rather,  we  intend  to  develop  in  this  paper 
the  properties  of  the  extreme  spread  more  extensively  than 
has  been  done  in  the  past,  and  thereby  contribute  to  an 
improved  understanding  of  the  statistical  characteristics  of 
the  probability  distribution  of  the  extreme  spread,  which  is 
required  in  any  first-class  or  overall  analysis  of  target 
accuracy  studies. 

In  our  introduction  of  the  subject,  we  point  out  that 

the  extreme  spread  is  a random  variable  which  follows  some 

kind  of  statistical  or  probability  distribution.  Indeed,  the 

amount  of  random  variation  from  one  group  of  shots  to 

another  depends  markedly  on  the  sample  size,  or  the  number 

of  shots  in  a group,  and  the  underlying  unknown,  population 

round- to -round  standard  deviation,  which  we  will  call 

o.  The  population  standard  deviation,  c,  is  a one-directional 

or  "linear"  quantity,  say  for  the  x or  horizontal  direction 

(as  well  as  the  y or  vertical  direction) , and  for  a very 

large  number  of  shots  it  may  be  found  as  the  square  root  of 

the  sum  of  squares  of  deviations  in  the  x-direction  from 

the  mean  divided  by  the  number  of  rounds.  In  rifle  firing, 

and  in  many  other  types  of  weapon  studies,  the  population 

standard  deviations  in  the  two  directions  are  equal  or  very 

nearly  so.  Hence,  it  may  be  assumed  in  >ur  following 

analysis  that  o « a ■ a.  The  exact  theoretical  probability 
xy 

distribution  of  the  extreme  spread,  or  bivariate  range,  as 
it  is  often  called,  has  not  been  determined  as  of  this  date, 
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although  many  o£  the  key  properties  of  the  distribution  are 
fairly  accurately  known  from  previous  studies.  In  the 
following,  we  report  on  the  results  of  a Monte  Carlo  type  of 
computer  simulation,  along  with  the  necessary  statistical 
analyses,  to  find  approximate  statistical  distributions  which 
for  all  the  practical  purposes  result  in  the  degree  of 
accuracy  needed  to  round  out  sufficiently  our  understanding 
of  properties  of  the  distribution  of  the  extreme  spread,  at 
least  for  the  very  important  practical  cases  involving  small 
sample  sizes  of  predominant  interest. 


Our  acknowledgements  must  go  to  Mr.  Philip  G.  Rust, 
retired  industrialist  of  the  Winnstead  Plantation, 

Thomasville,  Georgia,  for  his  great  interest  in  critical 
analyses  of  accuracy  firings  of  rifle;,  which  provided  much 
of  the  motivation  for  this  investigation,  as  well  as  for  the 
booklet  by  Grubbs  [1964],  which  are  of  importance  to  ballistic 
analyses  generally. 


II.  SOME  ANALYTICAL  PRELIMINARIES 


Consider  a random  sample  Xi Xn  (X^  = (x^  y^) 


from  a bivariate  normal  distribution  with  probability 
density  function  (p.d.f.)  given  by 

1 (x2+y2)/2o2 . 


f(x,  y) 


2irc2 
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The  extreme  spread  (ES) , or  bivariate  range,  is  defined  as 
ES  * max  |X. -X. |.  This  ES  is  of  course  a random  variable, 

• • 1 I 


as  previously  pointed  out,  and  we  seek  its  probability 
distribution,  realizing  that  it  will  be  dependent  upon  the 
sample  size  n.  For  the  case  n»2,  for  example,  we  have 


E3  « |Xi-X2j  * /(xi-X2)2+(yi-y2)2  a 2*o,  where  the  random 


variate  chi  has  two  degrees  of  freedom.  In  this  case,  the 
mean  value  of  the  extreme  spread,  E(ES)  **  1.77245o,  and  the 
variance  V(ES)  *(.9265o)2.  This  result  does  not  have  a 
direct  extension  to  higher  sample  sizes,  (n>2)  however,  and 
the  distribution  of  the  extreme  spread  ES  has  not  been 
determined  analytically. 


Some  earlier  work  of  Wilks  and  Grubbs  in  the  last 
reference  [1964]  have  led  to  Monte  Carlo  estimates  of  the 
first  four  moments  of  the  three-dimensional  or  trivariate 
range,  the  trivariate  midrange,  the  extreme  spread  or 
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bivariate  range,  and  the  bivariate  midrange.  Caeoullos  and 
DeCicco  [1967]  have  investigated  two  approximations  of  the 
ES  distribution  based  on  the  Monte  Carlo  moment  data  of  Wilks 
and  Grubbs.  In  this  paper  we  present  improved  Monte  Carlo 
moment  estimates  (by  virtue  of  greatly  increased  sample 
size),  as  well  as  quantile  estimates,  or  percentage  points, 
which  were  not  available  previously  and  propose  some 
approximate  distributions,  which  will  suffice  for  most 
analytical  studies  in  practice.  The  large-sample  moments 
tabulated  below  include  the  mean  u»  standard  deviation  o, 
the  skewness  measure  03,  and  "kurtosis"  or  peakedness  measure 
04.  (See  any  standard  textbook  on  statistics  for  further 

definitions  and  formulas.) 


III.  MONTE  CARLO  RESULTS 


The  moments  (Table  1)  and  various  quantiles  of  interest 
(Table  2)  are  based  upon  104  Monte  Carlo  samples  of  the 
random  variable  ES  for  each  value  of  sample  size,  n.  The  value 
of  n specifies  the  number  of  points  sampled  from  the  circular 
normal  distribution  to  determine  a single  value  of  ES. 

The  quantiles  of  Table  2 are,  of  course,  subject  to 
standard  jefinitions  and  interpretations.  For  example,  for  a 
samp*  size  of  n a 7 and  known  population  standard  deviation 
0,  the  lower  1*  point  is  1.842o  and  therefore  in  random  sampling 
from  a bivariate  normal  population  with  standard  deviation  o 
we  would  e-,ect  that  only  II  of  the  extreme  spreads  for  a 
sample  of  size  seven  would  fall  below  1.842o.  Similarly,  for 
the  991  point,  P ggQQ>  (or  upper  II  significance  level)  we 

would  expect  only  II  of  the  extreme  spreads  for  sample  size  n = 10 
to  exceed  S.75o.  The  mean  value  of  the  extreme  spread  for  a 
sample  of  size  10  from  Table  1 is  3.813o. 
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IV.  APPROXIMATE  PROBABILITY  DENSITIES 

Candidates  for  approximating  the  probability 
distribution  of  the  extreme  spread  were  chosen  from  three 
well  known  families  of  distributions,  specifically,  the 
chi  distribution,  the  lognormal  distribution,  and  the  two- 
parameter  Weibull  distribution.  Some  particulars  of  our 
findings  are  detailed  in  the  following  paragraphs. 

Chi  Distribution  Approximation 

In  considering  the  Chi  distribution,  we  made  use  of  the 
fact  that  2mw2/v  is  approximately  distributed  as  x2  with 


2m2 


degrees  of  freedom  (see,  for  example,  Grubbs  et  al 


[1966]),  where  wfl  is  the  sample  range  for  random  samples  of 

size  n from  a univariate  normal  population.  In  our 

notation,  m * E(w2)  and  v » V(w2).  We  used  this  same  type 

of  approximation  for  the  extreme  spread  ES  and  interpolated 
linearly  to  evaluate  chi-square  for  fractional  degrees  of 
freedom.  Using  the  sample  moments  of  ES  to  estimate  m and 
v,  this  family  provided  a rather  good  fit  to  the  sample 
auantiles  over  the  entire  range  of  n considered,  although  the 
dimension  of  the  sample  space  has  increased  from  one  to  two. 

The  results  are  summarized  in  Table  3 where  the  italicized 
value  is  the  fitted  value  juxtaposed  to  the  Monte  Carlo 
quantile  estimate. 

Lognormal  Distribution  Approximation 

The  lognormal  distribution  provides  an  excellent  fit  for 
large  values  of  n.  Following  recommendations  of  Aitchison 
and  Brown  [1966],  we  used  the  method  of  quantiles  (specific- 
ally, the  10th  and  90th  percentiles)  for  purposes  of  estimating 
v and  a of  the  associated  normal  distribution.  Possibly,  a 
different  choice  of  quantiles  might  lead  to  a better  fit  for 
small  n;  however,  a summary  of  results  for  the  larger  values 
of  n is  included  as  Table  4 with  the  same  format  as  Table  3. 

Weibull  Distribution  Approximation 

In  fitting  a two  parameter  Weibull  distribution, 

F(x)  ■ l-e*^0)  , to  our  data,  we  were  precluded  from 
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obtaining  maximum  likelihood  estimates  since  ve  have  at 
our  disposal  only  the  Monte  Carlo  moments  of  tne  distribution 
and  not  the  individual  sample  values.  In  lieu  of  these,  we 
used  a moment  estimate  suggested  by  Cohen  [1965],  but  with 
somewhat  less  than  satisfying  results.  In  fact,  both  the 
chi  approximation  and  the  lognormal  seem  preferable  to  the 
two  parameter  Weibull,  although  a three  parameter  Weibull 
with  the  introduction  of  a location  parameter  would 
probably  offer  some  improvement.  We  felt  that  this 
exceeded  our  charter  of  consideration  of  a few  commonly 
encountered  distributions  however,  and  postponed  any 
further  inquiry  for  a later  date,  especially  since  the  chi 
and  lognormal  distributions  gave  'ery  satisfactory  results. 


Data  Sot;:-  ry  and  T 1 air  le  of  Fitting  Procedure 

Th  chi  vaiia.f  is  less  dercriptive  in  terms  of  absolute 
difference  betwe.  . *'ne  Monte  Carlo  and  fitted  value  as  we 
go  further  out  in  tne  tails  of  the  distribution  and  the 
parameter  n increases.  As  a matter  of  fact,  the  upper  tail 
(perhaps  of  most  interest)  is  described  somewhat  better  than 
the  lower  tail,  although  the  percentage  error  between  the 
Monte  Carloed  and  fitted  value  rarely  exceeds  41  and  then 
only  in  the  extreme  percentiles  of  the  lower  tail.  It  is 
also  worthy  of  note  that  for  hypothesis  testing  the  region 
of  rejection  will  be  slightly  larger  than  that  indicated  for 
the  fitted  x variate. 

The  lognormal  variate  as  previously  stated  provides  a 
good  fit  for  the  larger  values  of  n(15  <_  n <,  34),  with  a 
percentage  error  in  excess  of  21  occurring  only  in  the  most 
extreme  percentiles. 

For  practical  situations  either  fit  is  adequate, 
the  chi  fit  being  more  versatile  over  the  range  of  n 
considered  and  the  lognormal  variate  offering  a closer 
approximation  over  a restricted  range  of  the  sample  size, 
n. 


Suppose  we  take  a random  sample  x^  . . . » x of  size 

n from  an  univariate  normal  distribution  and  determine  the 

sample  range  wn  » maix  |x. -Xj|.  Since  • w2  is 

1 » 3 


2m2 

approximately  distributed  as  x2  ( — M where  mj  « E(w2)  an<* 

vi 
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Vi  * V(w2),  we  wanted  to  see  if  the  extreme  spread,  ES, 

which  acts  like  a bivariate  sample  range  did  not  also 
closely  follow  a chi  distribution  (or,  equivalently,  ES2 

follow  a x2  distribution). 

To  -termine,  for  example,  the  95th  percentile  P 95 

of  ES2  corresponding  to  some  value  of  n we  must  satisfy 
the  relation 

Pr(ES2  < P>95}  - .95  * PrCES  < /F^T) 
or  equivalently, 

Pr{^B  ES2  < ~ P>95)  * .95 

To  test  our  approximation,  we  assume 

^ es2  . x2  (2si) , 

where  m and  v are  the  mean  and  variance  of  ES2,  so  that  we 
may  substitute  to  obtain 

Pr<x2(^>  < 22  P>95)  * .95 

and  nterpolate  in  the  chi  square  table  to  determine  P Q5. 

Finally,  multiplication  of  this  quantity  by  ^yields  P >95  for 
ES2,.  and  /P  g5  is  the  95%  point  for  ES. 

It  is  easy  to  show  that  the  mean,  m,  and  variance  v,  for 

the  extreme  spread  squared,  i.  e.  ES2 , may  be  expressed  in 
terms  of  the  moments  of  the  extreme  spread  ES  as  follows: 

m * a2  + v2 , 

and  v ■ 0404  + 4o3G3y  + 4c2v2  - o4» 

where  the  y,  0,  03  and  04  are  moments  of  ES. 
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Hence,  referring  to  Table  1 for  n = 10,  we  have 
m * (0 . 74S) 2 + C3.813)2  *=  15.094, 
v«  3.288(0. 7v5)  *♦  + 4 (0 . 388)  (0 . 745)  3 (3 . 813) 


I 

si. 

I 

e 

er 

I 


II 


+ 4(0.745)2  (3.813)2  - (0.745)*4  = 35.430, 
and  we  want  to  determine  P -,  where 

Pr(x2 (12 . 86)  <0.852  P g5>  » .95 
Interpolating  in  the  chi  square  table  for  12.86  d,  f.  yields 
0.852P  g5  « 22.20,  or  P g5  = 26.06  for  ES2 , and  /F^T  = 5.11 

which  is  the  95%  point  for  ES,  corresponding  to  the  entry 
in  Table  2. 

Example.  Compare  the  relative  precision  of  the  extreme 
spread  ES  and  radial  standard  deviation  USD  for  15  rounds. 


a 

For  n=*15  rounds  and  from  Table  1,  o = ES/4.190  gives  an 
unbiased  estimate  of  o,  and  quantity  .694/4.190  = .166 
is  the  relative  precision  for  the  extreme  spread.  In  a 
like  manner,  the  precision  of  the  RSD  is  found  from  Grubbs 
(1964)  Table  4 for  15  rounds  to  be  .1817/1.354  * .134. 

Therefore  .166  vs.  134  indicates  that  the  RSD  is  slightly 
more  precise  than  the  ES.  (For  the  relative  precision,  we 
compare  standard  errors  for  unbiased  estimates.) 

To  use  Table  2,  suppose  from  previous  firings  we  established 
that  a * 3 inches.  Then  for  n=15  rounds,  the  chance  that 
the  extreme  spread,  ES,  exceeds  5.396o  * (5.396)  (3)  = 16.19  inches 
is  .05. 


\ 


fe- 


ll 


' TV'  — ^ 


MOMENT  CONSTANTS  OF  THE  EXTREME  SPREAD 
TABLE  1 


n 

VES 

aES 

«3 

«4 

2 

1.766 

0.932 

0.632 

3.294 

3 

2.406 

0.887 

0.451 

3.143 

4 

2.787 

0.856 

0.393 

3.163 

5 

3.066 

0.828 

0.390 

3.171 

6 

3.277 

0.806 

0.374 

3.194 

7 

3.443 

0.783 

0.373 

3.177 

8 

3.582 

0.771 

0.392 

3.231 

9 

3.710 

0.754 

0.382 

3.215 

10 

3.813 

0.745 

0.388 

3.288 

15 

4.190 

0.694 

0.395 

3.255 

20 

4.452 

0.668 

0.400 

3.240 

25 

4.639 

0.650 

0.439 

3.307 

28 

4.734 

0.642 

0.426 

3.357 

30 

4.788 

0.635 

0.463 

3.441 

31 

4.822 

0.631 

0.434 

3.321 

34 

4.891 

0.623 

0.422 

3.318 

Note: 

The  numbers  in  the 

second 

column  are 

E(ES)/o,  and 

of  the  third  column  are  SD(ES)/c. 
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PERCENTAGE  POINTS  OP  THE  EXTREME  SPREAD 


TABLE  2 


n 

P.0005 

P.0010 

P.0050 

3 

0.339 

0.383 

0.578 

4 

0.653 

0.710 

0.946 

S 

0.885 

0.983 

1.260 

6 

1.137 

1.227 

1.491 

7 

1.348 

1.452 

1.710 

8 

1.525 

1.608 

1.863 

9 

1.607 

1.709 

2.030 

10 

1.798 

1.884 

2.140 

15 

2.295 

2.372 

2.656 

20 

2.630 

2.721 

2.972 

25 

2.894 

2.965 

3.220 

28 

2.952 

3.044 

3.312 

30 

3.084 

3.170 

3.402 

31 

3.149 

3.216 

3.429 

34 

3.216 

3.297 

3.517 

P.0100  P.02S0  P . 05 


0.687 

0.882 

1,.  066 

1 

1.076 

1.283 

1.477 

1 

1.400 

1.611 

1.801 

2 

1.636 

1.853 

2.043 

2 

1.842 

2.043 

2.243 

2 

1.998 

2.208 

2.403 

2 

2.167 

2.373 

2.563 

2 

2.277 

2.482 

2.669 

2 

2.772 

2.963 

3.129 

3 

3.095 

3.276 

3.438 

3 

3.329 

3.504 

3.652 

3 

3.424 

3.605 

3.759 

3 

3.511 

3.678 

3.834 

4 

3.541 

3.712 

3.868 

4 

3.630 

3,797 

3.946 

4 

CHI  APPROXIMATION 
TABLE  3 


P.0100 

P.0250 

0.687 

0.882 

0.676 

0.  868 

1.076 

1.283 

1.038 

2.  254 

1.400 

1.611 

l.  32S 

1.550 

1.636 

1.853 

1.560 

1.787 

1.842 

2.043 

l.  756 

1.984 

1.998 

2.208 

1.904 

2.  235 

2.167 

2.373 

2.060 

2.  255 

2.277 

2.482 

2.176 

2.404 

2.772 

2.963 

2.642 

2.864 

3.095 

3.276 

2.  953 

3.  Z73 

3.329 

3.504 

3. 171 

3.353 

3.424 

3.605 

3.284 

3.493 

CHI  APPROXIMATION 
TABLE  3 (CONTINUED) 


feS 

E&St 

|N 


fit' 


I 


I 


m 


I 


P.9Q00  P.9500  P.97S0 


P . 9900  P .9950  P .9990 


3 

3.588 

3.984 

4.318 

4.746 

5.002 

5.595 

3.583 

3.975 

4.306 

4.  717 

4.997 

5.576 

4 

3.916 

4.285 

4.602 

5.010 

5.290 

5.938 

3.923 

4.289 

4.597 

4.981 

5.  233 

5.  786 

5 

4.156 

4.519 

4.832 

5.207 

5.461 

6.057 

4. 164 

4.516 

4.813 

5.  177 

5.  424 

5 . ft3£ 

6 

4.336 

4.670 

4.973 

5.361 

5.655 

6.221 

4.352 

4.677 

4.966 

5.318 

5.55ft 

6.052 

7 

4.480 

4.805 

5.110 

5.471 

5.728 

6.245 

4.483 

4.799 

5.083 

5.417 

5.  647 

6.  227 

8 

4.595 

4.937 

5.227 

5.582 

5.848 

6.379 

4.608 

4.921 

5.  286 

5.521 

5.  752 

6.222 

9 

4.702 

5.029 

5.308 

5.672 

5.930 

6.398 

4.704 

5.012 

5.280 

5.  5ftft 

5.820 

6.277 

10 

4.786 

5.118 

5.409 

5.750 

6.004 

6.552 

4.  799 

5.105 

5.363 

5.680 

5.893 

6.343 

15 

5.101 

5.396 

5.668 

6.000 

6.235 

6.727 

5. 109 

5.387 

5.627 

5.915 

6.225 

6.524 

20 

5.336 

5.630 

5.880 

6.205 

6.436 

6.890 

5.333 

5.598 

5.833 

6.  265 

6 . 2i»3 

6.688 

25 

5.494 

5.790 

6.049 

6.364 

6.578 

7.012 

5.499 

5.  758 

5.980 

6.246 

6.42ft 

6.809 

2ft 

5.575 

5.860 

6.113 

6.453 

6.664 

7.138 

S.  582 

5.835 

6.056 

6.317 

6.497 

6.875 

30 

5.619 

5.898 

6.170 

6.476 

6.711 

7.205 

5.626 

5.877 

6. 100 

6.  35ft 

6.533 

6.905 

31 

5.651 

5.927 

6.180 

6.503 

6.719 

7.146 

5.654 

5.902 

6. 122 

5.  376 

6.553 

6.920 

34 

5.706 

5.979 

6.224 

6.523 

6.731 

7.218 

5.  715 

5.95ft 

5.  273 

5.424 

6.  5ft8 

6.958 

i 

j 

i 

i 
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LOG  NORMAL  APPROXIMATION 
TABLE  4 


n 

P.0005 

P.0010 

P.0050 

P.0100 

P.0250 

P.0S00 

P.1000 

15 

2.295 

2.397 

2.372 

2.478 

2.656 
2.  697 

2.772 
2.  an 

2.963 

2.986 

3.129 
3. 146 

3.340 
3.  348 

20 

2.630 

2.679 

2.721 
2.  761 

2.972 

2.984 

3.095 

3.098 

3.276 

3.274 

3.438 

3.433 

3.626 

3.826 

25 

2.894 

2.907 

2.965 

2.989 

3.220 
3.  222 

3.329 
3.  325 

3.504 

3.499 

3.652 

3.855 

3.845 

3.845 

28 

2.9S2 

3.020 

3.044 
3. 102 

3.312 
3.  323 

3.424 

3.438 

3.605 

3.609 

3.759 

3.785 

3.953 

3.953 

30 

3.084 

3.088 

3.170 
3. 170 

3.402 

3.391 

3.511 
3.  584 

3.678 
3.  676 

3.834 

3.830 

4.017 

4.017 

31 

3.149 

3.127 

3.216 

3.209 

3.429 

3.43C 

3.541 

3.542 

3.712 
3.  724 

3.868 

3.869 

4.055 

4.855 

34 

3.216 

3.202 

3.297 

3.284 

3.517 

3.S04 

3.630 

3.617 

3.797 
3.  788 

3.946 

3.942 

4.127 

4.127 
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